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ABSTRACT 

We present a joint weak lensing and X-ray analysis of 4 deg 2 from the CFHTLS and 
XMM-LSS surveys. Our weak lensing analysis is the first analysis of a real survey using 
shapelets, a new generation weak lensing analysis method. We create projected mass 
maps of the images, and extract 6 weak-lensing-detected clusters of galaxies. We show 
that their counts can be used to constrain the power spectrum normalisation erg = 
0.92+ 

oJo f° r ^™ — 0.24. We show that despite the large scatter generally observed 
in the M-T relation derived from lensing masses, tight constraints on both its slope 
£NJ ■ and normalisation M* can be obtained with a moderate number of so urces provided 

I> ' that t he covered mass range is large enough. Adding clusters from iBardeau et al.l 

CO \ (|2007l ) to our sample, we measure M* = 2.71 + °gi 1O 14 /i _1 M . Although they are 

dominated by shot noise and sample variance, our measurements are consistent with 
currently favoured values, and set the stage for future surveys. We thus investigate the 
dependence of those estimates on survey size, depth, and integration time, for joint 
' weak lensing and X-ray surveys. We show that deep surveys should be dedicated to 

the study of the physics of clusters and groups of galaxies. For a given exposure time, 
wide surveys provide a larger number of detected clusters and are therefore preferred 
for the measurement of cosmological parameters such as as and M*. We show that 
a wide survey of a few hundred square degrees is needed to improve upon current 
measurements of these parameters. More ambitious surveys covering 7000 deg 2 will 
provide the 1% accuracy in the estimation of the power spectrum and the M-T relation 
normalisations. 
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1 INTRODUCTION 



* Based on observations obtained with MegaPrime/MegaCam, 
a joint project of CFHT and CEA/DAPNIA, at the Canada- 
France-Hawaii Telescope (CFHT) which is operated by the Na- 
tional Research Council (NRC) of Canada, the Institut National 
des Sciences de l'Univers of the Centre National de la Recherche 
Scicntifiquc (CNRS) of France, and the University of Hawaii. This 
work is based in part on data products produced at TERAPIX 
and the Canadian Astronomy Data Centre as part of the Canada- 
France-Hawaii Telescope Legacy Survey, a collaborative project of 
NRC and CNRS. It makes use of photometric redshifts produced 
jointly by TERAPIX and VVDS teams, 
f E-mail address : joel.berge@jpl.nasa.gov 



In the currently-favoured hierarchical model of structure 
formation, clusters of galaxies have formed from the col- 
lapse of gravitational potential wells (e.g. IPeeblesl Il980l : 
|Padmanabhanlll993l : iLacev fc Cold 1 19931 : lLokadl200ll )~and 
are powerful cosmological probes. For instance, since they 
are sensitive to the expansion history of th e Universe, their 
abun d ance and spatial distribu ti on (e.g. 



120051 : 
| 2006h 

and their mass function (e.g. iLokas. Bode fc Hoffman||2004l ) 



19961: Wane & Steinhardtl 


19981: IHorellou & Berg 


Nunes. da Silva & Aghanim 


20061: iManera & Mot 
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depend on cosmological parameters, such as the dark en- 
ergy equation of state parameter w (e.g. Basilakos 2003; 
M aor fc Lahavll2005l : lBasilakos fc Voglidl2007|). or the power 
spect rum normalisation as (e.g. lSeliakll2002l ; |Pierpaoli et al.l 
120031 ). Several observational methods now permit the use of 
clusters of galaxies as cosmological probes, such as X-ray 
observations, weak gravitational lensing and the Sunyaev- 
Zeldovich effect. 

Due to improvements in telescopes and tech- 
niques, X-ray studies are able to constrain clus- 
ter physics and the mass scaling relation with ever 
greater precision. For i nstance, the self-similarity 
for clusters of galaxies jEke. Navarro fc Frenkl 1 19981 ; 
lArnaud. Aghanim fc Neumannlhood ) has been observation- 
ally verified. Nevertheless, adiabatic simulations still pre- 
dicts a mass-temperatur e relation with double the observed 
norm a lisation (e.g. Ncvalaincn, Markcvitch fc Formanl 



2000; iFinoguenov. Reiprich fc Bohringeij l200lf ) and the 
self-similarity assumption could break down at low 
temperatures (T ^ 3 keV). Thus, a steepening of the 
mass-temperature (M-T) relation is expected if galaxy 
groups underwent a preheati ng by supernovae, or a surge of 
entropy, in their early days (|Bialek. Evrard fc Mohrll200ll ; 
iMuanwong et al.ll2002l). Recent evidence for this steepenin g 
was found by e.g. iNevalainen. Markevitch fc Formanl |2000); 
IFinoguenov. Reiprich fc Bohringerl ( 20011 ). or Arnaud, 



Pointecouteau & Pra tt 2005 (APP05 hereafter), but coul d 
not be seen by e.g. lEttori. De Grandi fc Molendl d2002l); 



ICastillo-Morales fc Schindlerl l|2003h ; IVikhlinin et~aH 
(2006). Moreover, the M-T normalisation estimation is 
currently limited by systematics in measuring cluster 
masses from their X-ray profiles. This limitation can be 
lifted by using probes which are independent of the physical 
state of the cluster. 

Beyond galaxy cluster physics, the M-T relation is 
needed by X-ray experiments to estimate the power spec- 
trum normalisation as- Measuring this parameter has trig- 
gered much effort in several observational areas. For in- 
stance, CMB experiments fe.g. ISpergel et ahll2007l ) tend to 
a low value for as 0.8), weak lensing experiments tend 
to higher values 0.8), and X-ray observations provide 
intermediate values. 

Gravitational lensing does not depend on the underly- 
ing physics of clusters of galaxies or dark matter, but only on 
their potential wells, and on the Universe's geometry. Strong 
gravitational lensing has been used for galaxy clusters 
phys i cs (e.g. iMellier. Fort fc Kneiblll993l ; iKneib et al.|[l995l. 
19961 ; ISmail et al.l 1 19971 : iLuppino et al. Il999l; ISmith et ah! 
20051 ) and measurement of as (e.g. ISmith et al.l [2003). 
Weak gravitational lensing is more difficult to measure 
(for r eviews, see e.g. lMellierlll999l ; iBartelmann fc Schncidcr 
120011 ; iRefregierl l2003al ; iMunshi et ail I2006D). and has 
taken longer to be dete c ted jBacon. Refregier fc Elli si 



20001; Ivan Waerbeke et all |2000|; IWittman et all \20<M 



Rhodes. Refregier fc Grothll200ll ). Since then, particular at- 



tention has been given to cosmic shear, i.e. s t atistical cosmo- 
logica l weak lensing ( e.g. Bacon et al. 2003; Hcvmans et al.l 
20051; iMassev et all [20051. l2007al; iHoekstra et all 120061 ; 



Semboloni et alj|2006l ; ISchrabback et al.ll2007h . in attempts 
to measure w and as- It has also begun to be used a s 
a tracer of the cosmic web (e.g. IMassev et al. 2007b), 
and a way to detect and catalogue mass overdensi- 



ties (e.g. IWittman et all 120061. Gayazzi & Soucail 2007, 
GS07 hereafter, and Mivazaki et al.l [2007h ■ Beside the con- 
straints it can bring to cosmology, it can be used as 
a complement to X-ray analyses of clusters of galax- 
ies. Thanks to the physics-independent estimation of clus- 
ter masses, it appears as a unique method to calibrate 
the mass-temperature relation for clusters of galaxies 
e.g. Hjorth, Oukbir & van Kampenll998l;lHuterer fc White! 



I2OO2I ; iPedersen fc Dahldl2006l ; iBardeau et al.ll2007l ). It has 

been shown that the uncertainty in the normalisation of the 
mass-temperature relation is the largest source of erro r in as 
measurements inferred fro m X-ray cluster analyses (|SeliaJd 
120021 ; IPierpaoli et al.ll2003l ). An accurate mass-temperature 
relation, obtained from combined weak lensing and X-ray 
analyses, will thus provide new insights not only on the as 
discrepancy, but also on galaxy cluster physics. 



In this paper, we present the first joint analysis 
of weak gravitational lensing and X-ray wide-area sur- 
veys of a randomly-selected patch of sky. The weak lens- 
ing survey is derived from the CFHTLS, and the X-ray 
survey from th e XMM-LSS. The weak lensing analysis 
uses shapelets llRefregierl l2003bl ; iRefregier fc Baconl [2003; 
IMassev fc Refregierll2005l ). a new generation shear measure- 
ment technique, which has been shown to achieve a few per- 
cent ac curacy in shear measu rement from ground based tele- 
scopes (|Massev et al.ll2007d ). We analyse one square-degree 
of the CFHTLS Deep survey (the Dl field) and four con- 
tiguous square degrees of the CFHTLS Wide survey, which 
enclose the Dl field. We create convergence maps for this re- 
gion of the sky and give a catalogue of detections. We show 
how counting weak-lensing-selected clusters can provide an 
estimate of the power spectrum normalisation <T8- We then 
show how the combination of weak lensing and X-ray analy- 
ses of clusters provides an estimate of the mass-temperature 
relation normalisation T», independent of clusters physical 
state. Finally, we investigate the impact of a joint weak lens- 
ing and X-ray survey strategy on the accuracy of the as and 
T* measurement. We consider deep and wide weak lensing 
surveys, with the CFHTLS characteristics, combined with a 
blind X-ray survey of the same region of the sky. 



The organisation of the paper is as follows. Section [5] 
presents the surveys used in this paper, namely the CFHTLS 
and the XMM-LSS. The methods that we use are described 
in section 13 We show how we estimate the weak lensing ef- 
fect using shapelets, and how we generate convergence maps 
and detect clusters. We also briefly describe how the X-ray 
properties of clusters are obtained. Section [3] presents the 
convergence map we inferred from our weak lensing anal- 
ysis, and gives a catalogue of the galaxy clusters that we 
detect. We then give our estimates of the normalisation of 
the power spectrum and of the M-T relation. We then show 
in section[5]that combined blind surveys are necessary to get 
the best insights about those normalisations. The impact of 
survey strategy on those parameters estimations is discussed 
in section [6] We conclude in section [7] 
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2 DATA 

2.1 Weak lensing : CFHTLS 

The "Canada-France-Hawaii Telescope Legacy Survey" Q 
(CFHTLS), a joint France-Canada project, consists of three 
different surveys, namely the Very Wide Survey, the Wide 
Synoptic Survey (referred to as "Wide survey" hereafter), 
and the Deep Survey. Once complete, the Wide Survey will 
cover 170 deg 2 (divided into four distinct patches ranging 
from 49 deg 2 to 72 deg 2 ) in five filters (u*,g',r',i',z'), down 
to a magnitude i' « 24.5. Its main goal is the study of large 
scale structures by weak gravitational lensing and galaxy 
counts. The Deep Survey covers 4 different fields, each with 
an area of 1 deg 2 , in the same five filters, down to i' fa 28.5. 
It is primarily intended for Type la Supernovae studies but 
it is also useful for measuring large-scale structures. The 
CFHTLS images were obtained from observations with the 
MegaCam camera, made of a 36 CCD mos aic, of 2048 x 
4196 pixels each, with a 1 deg 2 field of view l|Boulade et al.l 
120031 ). 

In this paper, we present the weak gravitational lensing 
based mapping of 4 deg 2 of the Wide Survey (Wl patch), 
which include the 1 deg 2 field of the Deep Survey (Dl field), 
using both Wl and Dl images. The geometry of the fields 
that we use is shown in Figure [T] The data processing (as- 
trometry, photometric calibration, stacking of images) has 
been done by the CFHT community and TerapixQ. We use 
Wl images optimised for weak lensing : each field is the com- 
bination of 7 stacked images, each of 620 seconds exposure 
time. We use the T0003 release of the Dl field, consisting of 
275 stacked images, with a total 37.4 hour exposure time. 
The average seeing is 0.7 arcsec. We masked parts of the im- 
ages with saturated stars and/or too high a noise, by hand, 
so as not to bias our weak lensing results. This operation 
removes 10% of the original area covered by the data. We 
do not mask the ghosts created by spurious reflections on 
the telescope optics around saturated stars, but we eventu- 
ally remove the galaxies that they cover from our catalogues, 
since they are too noisy. The average galaxy density is 28 
arcmin -2 in the Dl image, and 13 arcmin -2 in the Wl im- 
ages. 



2.2 X-ray : XMM-LSS 

Designed to cover an area of several tens of square degrees 
up to a redshift z — 1, the XMM-LSS survey aims at de- 
tecting a significant fraction of clusters of galaxies, in or- 
der to constitute a sample suited to cosmological studies 
|Pierre et al.ll2004h . Its nominal exposure times are 10 ks, 
and have been raised up to 20 ks for the X MM Medium 
Deep Survey (XMDS : IChiappetti et al.ll2005l ). a 2 deg 2 re- 
gion included in the XMM-LSS, which covers the CFHTLS 
Dl field. In this paper, we use 4 deg 2 of the XMM-LSS which 
cover our 4 deg 2 CFHTLS data. The XMM pointings are 
shown on Figure [T] The raw X-ra y observations reduction is 
presented in iPacaud et ail i|2006l i. 



1 http://www.cfht.hawaii.edu/Science/CFHLS/ 

2 http://terapix.iap.fr 




Figure 1. Layout of the surveys. The red square is the CFHTLS 
Dl field. The four green squares are the four CFHTLS Wl fields 
used in this paper. Circles represent the XMM-LSS pointing s 
available in that region, prior to XMM A05 l lPacaud et 3.1.112007^ . 
Here, we only used those lying within the optical data. Those 
marked by a cross are strongly affected by flares and are being 
re-observed. 



3 METHOD 



3.1 



Weak lensing cluster detection 

Introduced inlRefrei der]|2003bl . iRefregier fc Baconll2003l and 



introduced m Kclrcgicr iUU3b, Kclrcgicr & Bacon 2U03 and 
iMassev fc RefregierlbOOa shapelets are a complete, orthog- 



onal, set of basis functions with which one can analytically 
decompose galaxy shapes. They can be seen as fundamen- 
tal shapes : a particular galaxy can be represented by a 
particular sum of shapelets basis functions Xn,m, weighted 
by coefficients f n , m - Their rich formalism provides an in- 
tuitive and analytic form for geometrical transformations 
(such as smear, shear, rotation) and for (de)convolution. 
Hence, they allow one to analytically describe the smear- 
ing of the Point Spread Function (PSF) and the shear of 
galaxies, properly correcting for the PSF. The shear estima- 
tion they provide has been shown to reach the needed ac- 
curacy for the CFHT LS specifications by the STEP project 
jMassev et al.ll2007cll . 

Our full pipeline will be described in an upcoming paper 
(Berge et al. in prep). We briefly summarise it here. Each 
sufficiently bright and non-saturated star is first decomposed 
into shapelets. A polynomial interpolation of each shapelet 
coefficient then provides a model of its spatial variations 
across the image. We are thus able to reconstruct the shape 
of the PSF at the position of each galaxy, the condition 
necessary for deconvolving it from the galaxies' shapes. Sev- 
eral stringent tests then validate our PSF model. In partic- 
ular, we require that the ellipticity distribution, and the two 
point correlation functions of the ellipticity, of the residuals 
between observed stars and their shapelet models are con- 
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sistent with zero. We also require that the cross-correlation 
between stars and galaxies ellipticity, when corrected from 
the PSF, is consistent with zero. The shape of galaxies is 
finally measured by decomposing them into shapelets, while 
decon volving from the PSF, as shown in lMassev fc Refregierl 
120051 . 

A shear estimator is created fro m the shapelet decom - 
positions of galaxies, as prescribed bv lMassev et al.|[2007dl : 



7 = 



/2,2 
Pi 



(1) 



where the shear susceptibility P 1 =< /o,o — /4,o > /v2 is 
fitted on the magnitude-size plane for galaxies. The coeffi- 
cients are complex numbers. The shear 7 of equation fT} is 
the complex notation for shears, 7 = 71 + 172. 

To increase the signal to noise ratio of our measure- 
ments, we give to each galaxy g a minimum variance weight 
w g — (crf s + cf> >g + cr?,,.) -1 , where o> i9 is the error on 
shape measurement for galaxy g, ap^ the error on the mea- 
surement of its shear susceptibility, and <7mt is the intrinsic 
ellipticity dispersion, set to <7i nt = 0.3. Slightly changing 
(Tint would be equivalent to giving more or less weight to 
our measurement errors, and would mostly affect the error 
bars in the shear measurement. The eventual peak detection 
would not be affected by such slight changes. This weighting 
scheme is equivalent to smoothly selecting the most useful 
galaxies for shear measurement. For instance, the faintest 
are down-weighted. It therefore provides us with effective 
densities of n e ff ~ 20 and 9 useful galaxies per square ar- 
cminute, in the Dl and Wl images, respectively. Then, a 
direct inversion in Fourier space of the pixelised shear map 
allows us to infer a con vergence (i.e. projected mass) map 
( Kais er fc Squired Il993h of the images. Structures in this 
mass map are extracted from the noise using a Gaussian 
filter. Figures [5] and [3] show the convergence maps that we 
inferred from our data. These are described in section 14.11 
While constructing a convergence map, we also create a 
signal-to-noise map, the signal-to-noise ratio being defined 
as 



v{x,y) = 



K{x,y) 
a K (x,y)' 



(2) 



where n(x,y) is the convergence at the (x,y) sky coordi- 
nates, and o K {x,y) its r.m.s error. The r.m.s error a K (x,y) 
is computed using Monte-Carlo simulations in which the in- 
put galaxies are positioned like the observed ones but with 
randomised shape orientations. 

Structures are then searched for in the filtered conver- 
gence map, a nd their astrometry prov ided, by the SExtrac- 
tor software (|Bertin fc Arnou"tall996r i. They are extracted 
according to their signal-to-noise peak, read from the signal- 
to-noise map. Hereafter, we define a 'significant structure' 
as a detection with a signal-to-noise ratio greater than 
2.5. Their mass is related to their integrated convergence 
through the lensing geometry, and can be estimated when 
their redshift and the redshift distribution of background 
galaxies are known. To account for the latter, we use the 
normalised distribution 



n(z) 



—?—(-)' 



exp 



(3) 



bv lBeniamin et all l|2007l ) (a,/3,z s ) = (0.836,3 425, 1.171) 
To ac count for n(z) in the Dl image, we fit lllbert et al.l 
(2006) 's photometric redshift distribution in the CFHTLS 
Dl field, and obta i n (a ,0,z a ) = (0.828,1.859,1.148). 
Ivan Waerbeke et al.l (|2006l ) have shown that errors in the 
n(z) fit are subdominant compared to Poisson noise and 
sample variance for the measurement of cosmological pa- 
rameters. We thus neglect them hereafter. 

We measure a cluster's virial mass by averaging its con- 
vergence in an aperture large enough that we can assume 
that the entire cluster is captured. The aperture corresponds 
to the region enclosed in the 2a level of the cluster's conver- 
gence map. This tech nique is similar to using a (^-statistic 
jFahlman et alJll994 ). with infinitely large annulus around 
the cluster, provided that the convergence in the entire field 
averages to 0. We verified this latter point, thus validating 
our choice. Note that because of the small number density 
of background sources, we cannot reliably fit a shear p rofile 
around clusters (see e.g. IPaulin-Henriksson et al.l (|2007T ) for 
an example of mass estimation using two profile fits around 
the galaxy cluster Abell 209). We then convert the virial 
mass into M2oo,c, the mass enclosed in the sphere of mean 
overdensity 200 times higher tha n the c ritical density, using 
the recipe from iHu fc Kravtsovl l|2003h . Hereafter, we will 
note M2oo,c more briefly M200. 

Weak lensing is affected by the entire mass distribu- 
tion along the line of sight. As a result, the weak lens- 
ing mass measurement of one cluster is biased by projec- 
tion effects. It has been shown, using different mass esti- 
mators, that large-scale structures in the line of sight, and 
near the ta rget cluster, introduce errors ranging from a 
few percent dReblinskv fc Bartelma"nnlll999l ; iHoekstrallioOll . 
l2003l:IClowe. De Lucia fc Kind2004l') to a few tens of percen t 
l|Metzler. White fc Lokenl l200ll ; Ide Putter fc White! I2005T I. 
In this paper, we assume that they produce a 20% error, 
added in quadrature to the shear measurement error. 

3.2 X-ray cluster detection and analysis 

The X-ray cluster d etection pipeline has been described in 
iPacaud et ail (2006). It takes account of the Poisson nature 
of the X-ray images, to extract and analyse clusters of galax- 
ies in a two-step procedure. Clu sters are first dete cted by 
a multi-resolution wavelet filter (jStarck et al.lfl99^ ). Then, 
each source is analysed using a maximum likelihood profile 
fitting procedure, and its X-ray properties a ssessed. Three 
classes of extended sour ces have been defined |Pacaud et alj 
l2006l : lPierre et al.ll2006h : (1) the CI class contains the high- 
est surface brightness sources, and is uncontaminated ; (2) 
the C2 class allows for 50 % contamination, and contains less 
bright extended sources ; (3) finally, the C3 class contains 
optically confirmed sources with apparent X-ray emission, 
that were not selected as CI or C2. In this paper, we only 
consider CI class detections, representative of the most mas- 
sive objects seen in the XMM-LSS. The redshift of detected 
clusters has been measured using spectroscopic observations 
from a number of tel escope and instrum ent combinations de- 
tailed in Table 2 of iPierre et alj (120061 ) . The ir temperature 
estimation is described in I Willis et all (|2005T ). 

IPacaud et al.l l|2007l ) have extracted and analysed 29 CI 



clusters from 5 deg of the XMM-LSS data (shown on Fig. 
where the parameters a, (3, z s are given for the Wide images [1]), that contain our 4 deg 2 optical data. Among other things, 
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Figure 2. Convergence map inferred from our weak lensing measurement of the Wl field. The square in the Wl map shows the boundaries 
of the Dl field (Fig. |3j. The map is smoothed by a 2.3' FWHM Gaussian. Dashed circles mark CI X-ray clusters. 



they have measured their luminosity and temperature. Here, 
we take into account their 16 clusters which are enclosed in 
the fields of our CFHTLS data, making use only of their 
temperature and redshift. They are listed in Table [1] Note 
that the cluster XLSSC053 is in the G12 XMM-LSS pointing 
(shown by the cross in the Dl fiel d, on Fig. [TJ , which was 
not used when lPacaud et all (120071 ) analysed the XMM-LSS 
observations. This pointing has been re-observed, and the 
X-ray characteristics of XLSSC053 are listed in Table [1] 



4 RESULTS 

In this section, we give the properties of clusters of galax- 
ies detected with our weak lensing pipeline. Counting the 
weak lensing detections allows us to constrain the matter 
power spectrum normalisation as. We then use the weak 
lensing mass of the detected groups to calibrate the mass- 
temperature relation for clusters of galaxies. 



4.1 Convergence maps and cluster catalogue 

Figure [2] shows the 4 deg 2 of the Wl field that we consid- 
ered. No significant overdensity (i.e. with v 2.5) has been 
detected. As we will quantitatively show in section 14.2.11 
this is consistent with the expected cluster counts for this 
survey. The black square in the image shows the position of 
the Dl field. Since there are around 20 galaxies per arcmin 2 , 
we expect more significant detections in this deep field (see 



section 14.2. 1 p . The upper panel of Fig. [3] shows the conver- 
gence map that we obtained from our weak lensing anal- 
ysis of the Dl field. Due to the varying level of noise in 
our map, which varies independently of re, two peaks with 
the same re value do not necessarily have the same signifi- 
cance. That results in the rejection of seemingly significant 
structures, such as the peak around (a, <5)=(36.75°, —4.75°). 
Significant structures are marked out by the white con- 
tours, which start at 2.2<r, with an increment of 0.5a. The 
lower panels show individual candidate clusters in more de- 
tails. In these, contours start at 2.2<r, with an increment 
of 0.2a. Even though we consider as significant a structure 
with v ^ 2.5, we plot the 2.2<r contours as a way to show 
the extension of our detect ions. D ashed circles mark X-ray 
cluste rs, and dotted circles rGSOTt 's weak lensing detections. 
GS07 measured the shear in the Dl field using the KSB 
method dKaiser. Squires fc BroadhurstHl995h . 

Table [T] summarises the measured characteristic s of th e 
clusters that we detect, together with all X-ray and iGSOTl 's 
detections in the region. The clusters that we detect through 
our shapelets weak lensing analysis are labelled with WLirf, 
where id runs from 00 to 07, and are sorted by decreasing 
significance. Their labels are listed in the first co lumn. Their 
official XMM names are given in column (2), and lGS07b IDs 
in the third column. X-ray clusters marked by a * are outside 
the Dl field. Columns (4) and (5) give their position. Col- 
umn (6) lists their spectroscopic redshifts, except for clusters 
WL03, for which a tomographic redshift is given, and WL04 
for which a photometric redshift is given (Aussel et al. in 
prep). The significances of the weak lensing detections are 
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Figure 3. Convergence map inferred from our weak lensing measurement of the Dl field. The convergence k is shown for the entire field 
on the upper panel. Lower panels show zooms on cluster candidates WL01, WL04, and the region surrounding WL03, WL05 and WL06. 
In the upper panel, contours levels start at 2.2<r, with an increment of 0.5cr. In the lower panels, they start at 2.2a, with an increment 
of 0, 2(7. Th e maps are smoothed by a 1.1' FWHM Gaussian. In the upper panel, dashed circles mark X-ray clusters, and dotted circles 
show lGS07l 's KSB weak lensing detections. Clusters detected by our shapelets weak lensing measurement are labelled WLid, and X-ray 
clusters not detected by weak lensing are labelled by their X MM name. All clusters are listed in Table[T] For clarity, false detections near 
edges are not shown. 
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Table 1. Clus ters catalogue. Besides the clusters that we detect through our shapelet weak lensing measurement, we also list clusters 
seen by IGSOTI and selected as CI clusters in the XMM-LSS. Weak lensing detections' significance is given for Dl and Wl, even if no 
detection appears in Wl. JV/20o(WL) is the cluster weak lensing mass. M2oo(X) is the mass estimated from X-ray profile, extrapolated 
from H500 to -R20O1 with respect to M500 masses given by Pacaud et al. (2007), and must be used with caution (see text). 



Weak lensing 


XLSSC 


GS07 


RA 


Dec 


z 


Significance 


M 200 (WL) 




A^oo(X) c 


Notes 


cluster ID 


number 


ID 


(deg) 


(deg) 




Dl/Wl 


(10 ld ft, i M Q ) 


(keV) 


nl Q 1 1 i r \ 

(10 iJ h l Mq) 




WL01 


013 


C103 


36.8497 


-4.5481 


0.31 


3.61 / - 


8.2+?i 

— 1.9 


1.0+°'l 

—0.1 


2.1 




WL02 






36.6589 


-4.7516 




3.09 / - 








WL03 




C104 


36.3628 


-4.1886 


0.32 a 


2.91 / - 


8.9+f'i 








WL04 


053 


C102 


36.1229 


-4.8341 


0.50 b 


2.90 / - 


10.3+|° 


3.4+3-1 

— l.U 


5.0 


XMM-LSS pointing 


















not observed in 


















1 3+0.1 




Pacaud et al. f2007) 


WL05 


041 


C114 


36.3723 


-4.2604 


0.14 


2.62 / - 




3.5 




WL06 


014 


- 


36.1389 


-4.2384 


0.26 


2.48 / - 


7 0+2.3 

' • -1.7 


1 q+0.2 


3.7 


just below detection threshold 




















in GS07's catalogue 


WL07 d 


022 


C107 


36.9167 


-4.8606 


0.29 


2.42 / - 


- 


L - '-0.1 


5.3 


near a mask 


- 
- 


025 


C105 
C110 


36.3375 
36.8167 


-4.6925 
-4.1269 


0.26 


-/- 
-/- 


- 
- 


2 0+°- 2 


6.5 


under a mask 




029 




■^R 01 79 


-4.2260 


1.05 


/ 
/ 




4 1+0-9 


13.9 


LUU J.llJiJ.1 I L.U.ollJ.1 L 




011 




36.5410 


-4.9680 


0.05 


-/- 




64+ 06 


1.0 






005 




36.7866 


-4.2995 


1.05 


-/- 




3.7+1 5 


16.5 


too high redshift 




006+ 




35.4382 


-3.7717 


0.43 


X/- 




4 8+°- 6 


30.4 


near an edge 




040t 




35.5232 


-4.5463 


0.32 


X/- 




1 6+ 11 
L -°-o.3 


6.8 






049t 




35.9892 


-4.5880 


0.49 


X/- 




9 9 +0.9 
z,z -0.5 


5.0 






018+ 




36.0079 


-5.0903 


0.32 


X/- 




2 0+°- 7 


8.0 






02lt 




36.2338 


-5.1340 


0.08 


X/- 




0.68_ Q Q2 


1.8 






ooit 




36.2378 


-3.8156 


0.61 


X/- 




9+0.4 
°- z -0.3 


14.3 






008 1 " 




36.3367 


-3.8014 


0.30 


X/- 




1 o+0.7 
■"•-^-O.a 


2.1 






002t 




36.3841 


-3.9198 


0.77 


X/- 




2 8+°- 8 


9.6 





a Tomographic redshift (GS07) 

b Photometric redshift (Aussel et al. in prep) 

c Rough estimates based on the isothermal assumption and extrapolated from M500 given by Pacaud ct al. (2007). 
d X-ray coordinates, 
t Outside Dl. 



listed in column (7), in Dl and Wl. A '-' means that the 
cluster is not detected ; a 'X' means that the cluster is out- 
side the Dl field. Columns (8) and (9) give their weak lens- 
ing mass Af2oo(WL) and X-ray temperature, respectively. 
Column (10) gi ves the mass e s timat e from the X-ray data, 
M 2 oo(X). As in iPacaud et~aH l|2007l ). these were evaluated 
under the assumption of an isothermal /3-model gas distribu- 
tion in hydrostatic equilibrium with the cluster's potential 
well. In this earlier work, the associated statistical errors 
were generally dominated by the temperature measurement 
uncertainty (8M/M w ST/T of order 10-25%). Here, the 
error on the emission profile can also become quite signifi- 
cant because we estimate the masses within -R200 instead of 
-R500 where t he X-ray emission start s to vanish. Moreover, it 
was shown bv lVikhli'nin et al.l l|2006l ) and lRasia et al.l (|2006l ) 
that the isothermal /3-model assumption leads to an under- 
estimation of the total mass, by up to 40% for low mass 
systems. For these reasons, we decided not to provide er- 
ror bars for our X-ray masses. Finally, column (11) gives 
some details about weak lensing detections, explaining for 
instance why we chose not to take them into account, or why 
we do not detect a cluster seen by another method. Among 
the rejection criteria are the proximity to an edge or to a 
masked region, the mass inversion procedure being sensitive 



to missing data and to edge effects. In that sense, and for 
clarity, the detections closest to edges have been removed 
from Fig. [2] and Fig. [3] 

Clusters WL01, WL04 and WL05 have unequ ivocal 
counterparts both in our X-ray catalogue and in lGS07r s KSB 
one. No significant CI X-ray sourc e has b een selected around 
WL02, and it remains invisible to lGS07l . Moreover, a visual 
inspection of the optical images does not show any galaxy 
overdensity around it. No significant CI X-ray source has 
been detect ed at t he position of WL03, even though it is 
also seen bv lGS07l . Cluster WL06 lies just below our detec- 
tion threshold {v = 2.48). Nevertheless, since GS07 detect 
it (though just below the detection threshold they use for 
their analyses) and since it is also detected by our X-ray 
analysis, we decided to list it, and to assess its weak lens- 
ing characteristics. The significance of cluster WL07 is even 
lower (y = 2.42). Since it is found close to the XLSSC022 
cluster (which coincides with lGS07l 's C107), we show its con- 
tours in Fig. |3] and list it in Table [T] However, it lies near 
an edge and a mask, so that its weak lensing characteristics 
are likely to be biased. We thus do not measure its mass, 
and will no t take it into account in what follows. Cluster 
XLSSC025 |GS07l 's C105) is under a mask, and cannot be 
detected by our pipeline. 
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In summary, out of our 7 shapelet weak lensing de- 
tections, 4 (WL01, WL04, WL05 and WL07) have a coun- 
terpart both in our X-ray catalogue and in the KSB weak 
lensing catalogue by IGS07I . even though we remove WL07 
from our subsequent analyses. One detection (WL02) ap- 
pears only in our catalo gue. O ne (WL06) has an X-ray coun- 
terpart, and appears in lGS07l 's map, but just below the de- 
tection threshold they use fo r their analysis. Finally, one 
(WL03) has a counterpart in GS07's catalogue, but is not 
selected as a CI X-ray cluster. This proves a good agreement 
between the three cluster detection methods used in those 
observations. X-ray clusters XLSSC005 and XLSSC029 are 
at too high a redshift to be detected with our surveys. Clus- 
ter XLSSC011 is too close and not massive enough to be 
detected by weak lensing, as will be shown in section 14.2.11 

Figure [4] compa res th e significance of our detections 
with that given by IGS07I . We make use of WL06, even 
though IGS071 did not use it, but gave its significance. 
Given our weak lensing measurement characteristics (n e e = 
19 arcmin -2 , a\ nt = 0.3) and theirs (n e g = 25.3 arcmin -2 , 
(Tint — 0.23), Eq. Q below allows us to compute the ex- 
pected proportionality factor between our detections' sig- 
nificance Kjhapdcts and theirs, ^gso7- We expect ^gso7 = 
1.47 ^shapciots- This relation is shown by the dashed line on 
Fig. [4] The significance of clusters in both catalogues scale 
as expected. One should note that this relation depends on 
the measurement characteristics for both methods in those 
particular experiments, and should not be used as a final 
comparison between KSB and shapelets. More comparison 
on real data will be needed in order to explore this issue. 

While X-ray masses listed in column (10) of Table 
[1] must be taken with caution, they can be compared to 
the weak lensing masses listed in column (8). Although 
one can notice an order of magnitude agreement between 
-M2oo(WL) and M2oo(X), masses estimated from X-ray data 
seem slightly underestimated. This is consistent with the 
previously mentioned bias arising from the isothermal beta- 
model parametrisation. 

While they do not provide strong statistics, our detec- 
tions can be used to estimate as, as shown below. Four de- 
tections have an X-ray counterpart and can thus be used 
to constrain the mass-temperature relation, provided that 
we add clusters from another catalogue. This is described 
below. 



4.2 Clusters number counts 

4-2.1 Weak lensing selection function 

A weak gravitational lensing selection function can be 
computed analytically (see e.g. Hamana. Takada fe Yoshidal 



|2004| ; iMarian fe Bernsteinll2006 ) from the signal-to-noise ra- 
tio of a halo in a given cosmology and weak lensing sur- 
vey parameters. We derive such a selection function, using 
an optimal match filter, in Berge, Amara & Refregier (in 
prep). In an observation characterised by a number density 
of background galaxies n g , an NFW halo of convergence k 
has signal-to-noise ratio : 



6 
5 

4 r 



<P 1 



° 3 



<^L04 / 



<3/VL03^/VL01 
<-Wl_'07 

MIL06 



Figure 4. Comparison between our detection significanc es and 
those of lGS07l , Symbols are data points. C110 is detected bv lGS07l . 
but remains invisible to our pipeline. The dashed line features the 
expected relation between clusters' significance in both analyses, 

V GS07 = 1-47 ^shapelets- 



and 



<j 7 is the r.m.s shear error per galaxy, 
we neglect projection effects and sample vari- 



d 2 x k 2 (x) 



(4) 



where 
where 

ance, which have been sho wn to have subdominant effects 
ijMarian fe Bernsteinll2006h . 

Our selection function is shown in Figure in the 
mass-redshift plane, for our Deep (n 9 = 20 arcmin" 2 , 
a 1 = 0.3, thick black) and our Wide (n g — 9 
arcmin -2 , a 7 = 0.4, red) surveys, in a cosmologi- 
cal model based on the three-year Wilki nson Microwav e 
Anisotropy Probe results (WMAP3; ISpergel et"al] 120071 ). 
(/i,fi m /i 2 ,fi b /i 2 ,a 8 ,^)=(0.73,0.127,0.0223,0.76,-l). We use 
the redshift distributions for background galaxies given by 
equation (|3}. 

Figure [5] shows, from bottom to top, the minimum de- 
tectable mass for a halo at a 2a, a 3a and a 4a detection 
threshold. The Deep and Wide selection functions have dif- 
ferent slopes, as a consequence of their different redshift 
distribution for background galaxies. Symbols represent the 
position, in the redshift-mass plane, of clusters listed in Ta- 
ble [T] We use the weak lensing mass M2oo(WL) for WL01, 
WL03, WL04, WL05 and WL06 (thick square symbols) and 
the X-ray mass M2oo(X) for other clusters (diamonds). Al- 
though we detect cluster XLSSC022 (WL07), we do not as- 
sess its gravitational mass, and thus show its X-ray mass 
in Fig. [5] C luster XLSSC025 should be detectable (and is 
detected bv lGSOTT ), but it is under a mask in our analysis. 
Triangle symbols (labeled ' in Table [TJ correspond to CI X- 
ray clusters in Wl that are outside the Dl region. It is clear 
from Figure [5] that they can not be detected by our weak 
lensing analysis of Wl. Only XLSSC006 should be seen, at 
the 2a level. However, its detection is plagued by its prox- 
imity to the edge of the image. 

Also shown on Fig. [S] are the X-ray selection functions, 
for 50% and 80% detection prob abilities (dashed lines) as 
estimated bv lPacaud et alj (|2007l ). 
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Figure 5. Weak lensing selection function for a survey like Dl 
(thick black ; cri nt = 0.3, n 9 =20 arcmin -1 , redshift distribution 
as Eq. J3}) and Wl (red ; <r; nt = 0.4, n g =9 arcmin" 1 , redshift 
distribution as Eq. (0), assuming a WMAP3 cosmology in each 
case. From bottom to top, lines correspond to 2<r, 3cr and 4<r 
significance. Dashed lines show the X-ray selection function, cor- 
respo nding to 50%, and 80% detection probability jPacaud et al.l 
120071 Fig. 18, lower and upper curves, respectively). Thick square 
symbols arc our detections in the Dl data, labeled by their ID ; 
they are not detectable in the Wl data. Diamonds are clusters de- 
tected either by |GS07l or by X-ray analysis, in Dl, that we do not 
detect for reasons listed in the text. Red triangles are CI X-ray 
clusters lying outside the Dl region. Except for the thick square 
symbols (for which we use the weak lensing mass M20o(WL)), 
we use the X-ray mass Af2oo(X). Except XLSSC006 (labeled for 
visibility as X006), they cannot be detected by a weak lensing 
experiment in the Wl data. XLSSC006 is not detected because 
of its proximity to an edge. 



Figure[5]shows an excellent agreement between the clus- 
ters characteristics and their predicted detectability by weak 
gravitational lensing. 



4.2.2 Number counts 

From equation Q, the expected number of haloes detected 
above a certain signi ficance can be computed using a Press- 
Schechter approa ch JPress fc Schechterlll974i ). For this pur- 
pose, we use the I Jenkins et al.l (|2001 ) mass function to es- 
timate the number of haloes that we can detect, as a func- 
tion of significance threshold. Curves on figure [6] show such 
count s for different ag and survey depths. iMivazaki et al.l 
(2002) already used this statistic to discriminate between 
halo profile models. It is used here to measure ag. 

Most of our detections are validat ed by corresponding 
objects either in the catalogue of lGS07l or our X-ray CI clus- 
ter catalogue. Nevertheless, despite its relative high signifi- 
cance, WL02 does not have such independent support. We 
indeed consider it as a false detection, and do not take it into 
account for cluster counts. We then estimate the number of 
false detections from Monte Carlo simulations. For this pur- 
pose, we input galaxies at the position of the actual ones, 



Figure 6. Cumulative cluster number density as a function of 
weak lensing detection significance, in the Dl data. The error bars 
include shot noise and sample variance. The dashed line shows the 
expected number counts in a WMAP3 cosmology, for the survey's 
characteristics. The solid line is our best fit, when varying erg 
(erg = 0.92). The dot-dashed line is the expected number counts 
for the Wide survey, with erg = 0.92. 



but randomise their shear, and look for detections with sig- 
nificance higher than 2.5a. The convergence maps that we 
infer from them show only false detections. We find that, in 
this particular experiment, we expect only one false detec- 
tion above the 2.5a level. This is thus consistent with re- 
moving WL02. For this counting purpose, we remove WL06 
and WL07 from our catalogue, since they do not reach the 
2.5a level. The symbols on Fig. [6] represent our cumulative 
counts, corrected from false detections. Their error bars in- 
clude the effects of shot noise and sample variance, computed 
from Hu & Kravtsov's (2003) analytic formula. We then fit 
the expected number counts to our data as a function of as, 
keeping all other parameters constant. In order to avoid co- 
variance between our data points in the cumulative counts 
depicted by Fig. [S] we performed the fit on the expected dif- 
ferential number counts dN/du(v). We find ag = 0.92^^ 
(at the 68.3% confidence limit), for Q m = 0.24. Despite large 
error bars, we can set interesting constraints thanks to the 
strong dependence of these counts on ag, as shown by the dif- 
ference between the solid and dashed curves in Figure|6] The 
dashed line shows the expected cumulative number counts 
for the Deep survey in a WMAP3-like universe (ag — 0.76). 
The solid line is our best fit (ag = 0.92). The dot-dashed 
line shows the expected number density on the Wide sur- 
vey, with ag = 0.92. 



4.3 Mass-temperature relation 

Under the virial equilibrium assumption, the mass and tem- 
perature of_a_cIustCT_are_j^lated by the following scaling 
relation (|Pierpaoli et al.ll2003l ) 
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3ardeau et al (2007) 
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Figure 7. Mass-temperature relation, normalised to z = 0, for 
our group sa mple (diamonds). To i mprove the statistics, we added 
clusters from IBardeau et al.ll20oH (triangles). We use X-ray tem- 
perature and weak lensing mass. The so lid line is our best fit 
(Eq. 0. The dashed and dotted line are IAPP051 M-T relation, 
when they consider all clusters or only those with T > 3.5 keV, 
respectively. 



M vir (T, 
W 5 h- 1 M ( 



© 



3/2 



-1/2 



,0 A (z) 



where M v i r is the virial mass, T is the virial temperature, 
T» is a normalisation factor, and E(z) 2 = Q m (l + z) 3 + 
Q.A + + 2) 2 • Ac (z) is the overdensity inside the virial 
radius, in units of th e critical density. We compu t e it u sing 
the fitting formula by I Weinberg fc Kamionkowskil (|2003h for 
Avir = A e / n m , which i s very similar to an earlier approx- 
imatio n by iNakamural |l996) given by iKitavama fe Sutol 
l|l996h for a universe with arbitrary Q m . 

A more general relation often used to fit observations 
makes use of a related normalisation factor M* and is given, 
at redshift z — 0, by 



(6) 



where M200 is the mass inside the sphere of mean overden- 
sity 200 times higher than the critical density and a = 3/2 in 
the hydrostatic equilibrium assumption (e.g. IAPP05I ). Here- 
after, to account for redshift evolution, we normalise all our 
temperatures to z — by dividing them by E(z) 2 ' 3 . 

Measuring erg from X-ray c ounts is affected b y 
the degeneracy Q m as oc rr 8 l|Pierpaoli et all [2003). 
iPierpaoli et~afl (|2003l ) have shown that the uncertainty in 
M, is the main concern in measurements of as from X-ray 
cluster observations alone. Such data is limited by the re- 
quirement t hat the cluster mas ses be inferred from the X- 
ray profiles. ISmith et al. (2003) have also shown that un- 
relaxed clusters, being hotter than relaxed clusters, pro- 
vide a supplementary bias to the as estimate. It is thus 
important to have a mass estimate independent of the hy- 
drostatic equilibrium assumption. Weak gravitational lens- 
ing gives such an estimate. Combined with X-ray tem- 



perature, it can be efficiently used to constraint the M- 
T relation, independently of t he clu st er physical state. 
iHiorth. Oukbir fc van Kampenl <|l99Sl ); iPedersen fe Dahld 



(2006); iBardeau et alj (|2007f ) have already used it to mea- 



sure the M-T relation normalisation. 

As described above, we have the weak lensing mass 
and X-ray temperature of only four groups. Hen ce, to in- 
crease our statistics, we add lBardeau et all (|2007l Vs clusters 
to our catalogue, pro viding us with 11 additional clusters. 
IBardeau et alj (|2007T ) estimated cluster masses by fitting an 
NFW model to their tangential shears. We should note here 



that since IBardeau et"alT(|2007l Vs mass estimation and ours 
are based on different techniques, our subsequent analysis of 
the mass - temperature relation could be slightly biased due 
to possible calibration differences. Figure [7] shows the re- 
lation between the temperature and the weak lensing mass 
A/200 for the combined catal ogues. Diamonds are our groups, 
labelled WLi rf, triangles are Bardeau et all (|2007l )'s clusters, 
labelled Aid IBardeau et al. (2007) proceeded to the weak 
lensing analyses of massive halo e s, the temp erature of which 
were obtained by IZhang et al.1 l|2007h and lOta fc Mitsudal 
(2004) ; particularly, they estimated their weak lensing mass 
and measured the scale relations for those clusters. They 
fitted their sample by varying both 01 and M„, and found 
a large slope, far from the hydrostatic equilibrium assump- 
tion, a = 4.6 ± 0.7. Doing the same analysis on the larger 
range in mass that the addition of both catalogues probes, 
from galaxy groups to galaxy clusters, we find : 



-3/2 



(•») 



M 2 , 



lO 14 /i- 1 M 



2 ?1 +o.79 ( T V 
2 ' 71 - - 61 U kevj 



(7) 



which is in good agreement with lAPP05l (w hether they use 
all clu sters or only the most massive ones) . | Bardeau et al.l 
J2007l) (when th ey fix a = 1.5). |Pedersen fc Dahld l|2006l ) or 
iHoekstral l|2007i ). The solid line on Fig. [7] is our best fit. The 
dashed line is the best fit from IAPP05I . when they consider 
all clusters in their catalogue. The dotted line is their best 
fit when they consider their most massive clusters (T > 
3.5 keV). 

Our result is consistent with self-similarity evolution for 
galaxy clusters down to low temperatures. It is also consis- 
tent with previous measurements which observe a steepening 
of the M-T relation at the low mass end, due to the expected 
self-similarity breaking for such masse s (e .g. Nevalainen et 
al. 2000; Finoguenov et al. 2001; |APP05D . Moreover, one 
must be aware that the galaxy groups we consider were de- 
tected just above our weak lensing selection function (Fig. 
0. Due to the expected scatter in the M-T relation, those 
groups can represent only the most massive ones with tem- 
perature ranging about 1 keV. Our group sample could thus 
bias our fit towards a flat slope for the M-T relation. The 
analysis of more low temperature groups will be needed to 
further explore this issue. 



5 DISCUSSION 

The power spectrum normalisation as has been measured 
with different probes, such as X-ray clusters of galaxies, 
CMB, and cosmic shear (i.e. statistics of weak gravitational 
lensing) . Some discrepancies have emerged between the pre- 
ferred value from those measurements. Recent CMB ob- 
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Figure 8. Domain allowed for by different measurements, in the 
T» - cr8(Sl m /0.24) 6 plane. The shaded region shows the con- 
straints given by our <rg and T* measurements. Thick lines are 
our best fits. The slanted black region correspon d to the ler con- 
strain ts on the Q^ 6 crg oc T^~ ' 8 relation from IPierpaoli et al.1 
12001 The vertical, light blue, shaded region show s the la error 
011 T* from lAPPOa CMB derived constraints of erg l lSpergel et al.l 
l2007t ) are shown by the horizontal dark blue shaded re gion. 
Cosmic shear as es timations from IHoekstra et al J |2006| and 
IBeniamin et al.ll2007| are marked by the red and green horizon- 
tal shaded regions. iHoekstra et alj l2006's constraints are typical 
of cosmic shear results. They are hig her than X-ray estima tions, 
marked by the intersection between IPierpaoli et al.l |2003| 's and 
lAPP05l 's allowed domains. 



made with CFHTLS Wide data (red). This highlights the 
discrepancy between X-ray an d weak lensing estimat es of 
erg mentioned above. However, IBeniamin et al.l |2007) give 
a lower estimate for as, consistent with X-ray measurements 
(green). This could be a sign that other cosmic shear anal- 
yses did not take some systematics into account, and hav e 
thus overestimated erg. According to lBeniamin et all l|2007h . 
previously published analyses made use of i nsufficiently- 
accura te galaxy photometric redshifts. Using lllbert et al.l 
(2006) 's redshifts yielded a l ower value of erg both for cosmic 
shear (|Beniamin et al.ll2007l ) and for our cluster count anal- 
ysis. We found a 5% decrease in our as esti mation when 
going from previous redshift distributions to lllbert et al.l 
(2006)'s ones. This i s less than the change reported by 
iBeni amm et al.1 (|2007l ). and our best fit still tends to favour 
a higher va l ue fo r erg, but is limited by low statistics. 
ISmith et al.1 1120031) have analysed the bias from unrelaxed 
clusters in erg measurement using lensing clusters and the 
M-T relation. They found that unrelaxed clusters are 30% 
hotter than relaxed clusters : using unrelaxed clusters can 
provide 20% overestimates of erg. This is enough to explain 
the large range of measured erg, from ~ 0.6 to ~ 1. They es- 
timated crg(f2 m /0.24) ' 6 = 0.86±0.23. Estimates from X-ray 
alone can also be affected by systematics, such as the mass 
estimate from X-ray profiles of clusters. For example, a slight 
decrease of T* would cause an increase of the X-ray estimate 
for as- A better insight into this will come from an accurate 
measurement of T*, preferably with mass estimation meth- 
ods independent of cluster physics. Large combined weak 
lensing and X-ray surveys will be needed to disentangle the 
situation. They will provide both independent constraints 
on erg, and insights on T*. 



servations favour a low erg and cosmic shear us ed to em- 
phasise a high value (see e.g. iRefregierl l2003al for a re- 
view). X-ray clusters provide intermediate measurements 
(see e.g. IPierpaoli et al.l 12003 for a review). The dominant 
discrepancy between cos mic shear an d X-ra y clusters has re- 
cently been reduced by I Jarvis et al.l ((2006), who measured 
erg « 0.81 f or f2 m = 0.26 when u sing cosmic shear alone, 
follow ed by IBeniamin et al.1 l|2007l ). who used lllbert et al.l 
(2006)'s improved galaxies photometri c redshifts, and mea- 
sured erg = 0.84 for Q, m = 0.24, and iFu et all (|2008T l who 
found a consistent value. To clarify these discrepancies, one 
needs to measure both the power spectrum and the M-T 
relation normalisations erg and M*, as we discuss here. 

The regions allowed for by different measurements on 
the T» - erg(n m /0.24) o ei plane are shown by Figure [8] The 
shaded region shows the constraints given by our erg and T, 
measurements. Our best fits are shown by the thick lines. 
The slanted black band on the fig ure is the 68.3% boun d on 
the f^ 6 crg oc T' - 8 relation from IPierpaoli et all (|2003l ) us- 
ing X- ray clust ers. Its intersection with the vertical light blue 
band |APP05l 's T, estimation) gives the current value for 
as favoured by X-ray cluster o bservations, erg » 0.7 7 ± 0.06 
for an Q m = 0.3 universe fe.g. IPierpaoli et al.ll2003T ). which 
corresponds to crg(f2 m /0.24) ' 6 ~ 0-88 ± 0-05 . This value 
is higher than that measured bv lSpergel et al.1 ([2007) from 
CMB analyses of WMAP3 (dark blu e), but lower than most 
cosmic shear analysis, like that of IHoekstra etHI (|2006T l 



6 PROSPECTS FOR FUTURE SURVEYS 

In the following, we investigate the impact of future com- 
bined blind weak lensing and X-ray surveys on the mea- 
surement precision of the power spectrum and the mass- 
temperature relation normalisations. We take the WMAP3 
l|Spergel et al.l 120071) cosmology as our fiducial model. We 
consider two different ground based survey strategies for 
our weak gravitational lensing analysis : deep and wide sur- 
veys similar to the CFHTLS Deep and Wide surveys. We 
use their observed weighted number density of useful back- 
ground galaxies to be n g = 20 arcmin -2 and 9 arcmin -2 , re- 
spectively, distributed according to Eq. Q. We also assume 
the intrinsic ellipticity and shape measurement error to be 
er 7 = 0.3 in both cases. Following the CFHTLS scheduling, 
we take for exposure times 40 hours per square degree for 
the deep survey and 1 hour per square degree for the wide 
survey. 

6.1 as measurements 

We first investigate the impact of future surveys on the 
erg measurement. Using the Press-Schechter approach de- 
scribed in section 14.2.21 we estimate the number of weak 
lensing detections with significance higher than 2.5, taking 
into account shot noise and sample variance. We assume 
that all clusters have a spherically symmetrical NFW pro- 
file. We thus neglect the effect of haloes' asphericity shown 
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Figure 9. Relative errors on as, from clusters counts in a weak 
gravitational lensing survey, as a function of survey size and in- 
tegration time. All other parameters are kept constants. We as- 
sume that 1 deg 2 of wide requires 1 hour of observation time, 
and 1 deg 2 of deep requires 40 hours of observation time. That 
is, the lower x-axis shows area as well as the wide survey ex- 
posure time ; the upper x-axis shows the deep survey exposure 
time. The thick solid line corresponds to a deep survey, and the 
thick dashed line to a wide survey. The flat lines show the current 
error measurement from cosmic shear statistics (dash-dot, Hoek- 
stra et al. 2006; Benjamin et al. 2007), and from X-ray clusters 
(dash-dot-dot, Pierpaoli et al. 2003; APP05). 

by Clowe et al. (2004) : triaxial haloes oriented along the 
line of sight appear more massive than triaxial haloes of the 
same mass, but perpendicular to the line of sight, and thus 
have a higher signal-to-noise ratio. Clowe et al. (2004) have 
shown that this approximation does not yield any difference 
in the mass measurement dispersion. Figure [5] shows the 
68.3% relative error on <T8 that can be reached by counting 
weak lensing detected clusters as a function of their signifi- 
cance, for a deep (thick solid line) and a wide (thick dashed 
line) surveys, as a function of survey's size and observing 
time. Because of the higher number density of clusters it 
allows one to detect, a deep survey provides errors 3 times 
lower than a wide survey of the same size. However, for a 
given exposure time, a wide survey provides errors 2.1 times 
lower than a deep one. That means that the gain due to 
the coverage (and detectable clusters number) increase is 
faster than the one due to depth increase. A larger cover- 
age is also advantageous in that it makes sample variance 
fall down rapidly. Moreover, a wide survey detects the most 
massive haloes, the physics of which is better understood. 
Consequently, in a survey strategy driven by exposure time, 
one should prefer a wide survey. The flat dashed-dot line 
shows the current constraints provided by cosmic s hear anal- 
yses (jHoekstra et al.ll200fjl : iBeniamin et al.ll2007n . The fiat 
dashed-dot-dot line shows the current c onstrain ts from the 
combination of X-ray M, measurem ent (|APP05l ) and X-ray 
cluster counts (|Pierpaoli et al.ll2003T l. 

Detecting and counting clusters on a 10 deg 2 deep sur- 
vey will be competitive with current cosmic shear measure- 
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Figure 10. Relative errors on T* from combined weak gravita- 
tional lensing and X-ray surveys, as a function of survey size and 
integration time. All other parameters are kept constants. The 
lower x-axis shows area as well as the wide survey exposure time; 
the upper x-axis shows the deep survey exposure time. Solid lines 
show the errors for a deep survey, assuming the fractional error 
on weak lensing mass measurement is 20% (thick black) and 30% 
(red). Dashed lines, with the same colour indexing, show the er- 
rors for a wide survey. The current error measurement from X-ray 
clusters is shown by the flat dashed-dot-dot line (Pierpaoli et al. 
2003; APP05). We made the same assumption about the relation 
between survey area and observation time as in Fig. [9] 

ments, whereas 20 deg 2 of coverage is needed to compete 
with current X-ray clusters measurements. Those figures 
transform as 100 deg 2 and 200 deg 2 for a wide survey. That 
is, to compete with current cosmic shear surveys, one needs 
400 hours of deep survey exposure, or 100 hours of wide sur- 
vey exposure. Double these times are required to compete 
with X-ray surveys. A wide survey, less demanding in expo- 
sure time than a deep one, should then be used. Counting 
clusters on the entire planned CFHTLS Wide Survey 170 
deg 2 will provide a 6% fractional error on the as measure- 
ment. Reaching the 1% fractional error will require a 7000 
deg 2 wide survey, or a 700 deg 2 deep survey. Future surveys 
(e.g. Pan-STARR93, LSSlQ, DUNEO) will be able to achieve 
such errors. 

6.2 T» measurements 

We now turn to the precision that can be reached on T* 
measurements by future joint surveys. What matters here 
is not cluster counts as a function of significance, but as a 
function of mass. As we see from figure [5] low mass clus- 
ters of galaxies cannot be seen through weak gravitational 
lensing since they do not create high enough signal-to-noise 
ratios. A deep survey captures lower mass clusters than a 

3 http:/ /panstarrs. ifa.hawaii.edu 

4 http://www.lsst.org 

5 http://www.dune-mission.net 
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wide survey, but deep and wide surveys give access to the 
same number of massive clusters. Therefore, a deep survey 
is naturally focused on the physics of galaxy groups (e.g. it 
can probe similarity breaking at the low mass end of the 
mass-temperature relation). A wide survey gives the same 
statistics on massive haloes, generally used to measure the 
mass-temperature relation normalisation : for this purpose, 
one should then choose a wide survey. To compare the merits 
of both deep and wide survey on the T„ estimation's preci- 
sion, we simulate M-T relations for both types of survey. We 
take a realistic scatter int o accoun t, a\ 0&t iat = 0.051 for the 
logarithmic M-T relation i|APP05h . We assume that masses 
are measured through weak gravitational lensing. We mea- 
sure T* and the error on its estimate, by assuming (1) that 
our cluster sample is complete, (2) that we only make use 
of those clusters detected in our blind survey, and (3) that 
we know the X-ray temperature of each of them. We also 
investigate the influence of the mass estimation fractional 
error 5M/M. 

Figure [10] shows the 68.3% error on T* that can be 
reached from a combined blind X-ray and deep (solid line) or 
wide (dashed line) weak gravitational lensing surveys, as a 
function of survey's size and integration time. Here again, as 
for the error on the power spectrum normalisation, a deep 
survey gives errors 2.3 times lower than a wide one with 
the same sky area coverage. On the other hand, a wide sur- 
vey gives errors 2.7 times lower than a deep one with the 
same exposure time. The dependence on area underlines the 
T„ estimation's reliance on the number of useable haloes 
for the M-T relation fitting. Figure [TD] also shows the sen- 
sitivity of the T t estimation to the mass measurement er- 
rors. The black lines assume SM/M — 0.2, and the red ones 
SM/M = 0.3, which are the current fractional errors on mass 
measurement from weak lensing. Going from SM/M = 0.3 
to SM/M — 0.2 allows one to reduce the error on T* by 
a factor of 1.3 (resp. 1.2) for a deep (resp. wide) survey of 
a given sky area. The flat solid line repr esents the current 
error on T* from X-ray clusters (|APP05l ). Assuming a 20% 
error measurement on weak lensing masses, one needs a 50 
(resp. 300) deg 2 weak lensing deep (resp. wide) survey to 
reach the current error. Reaching the 1% fractional error 
(for our fiducial model with T* = 1.9) will require a 2500 
deg 2 wide survey, or a 500 deg 2 deep survey. Weak lensing 
surveys like LSST or DUNE combined with X-ray surveys 
like eROSITA will be able to reach such a limit. 

In a survey strategy driven by exposure time, a wide 
survey of 2500 deg 2 (2500 hours) will be able to reach the 
1% accuracy both on erg and T» , at a much cheaper expense 
than a deep survey. Nevertheless, a deep survey will still be 
useful to probe high redshift regions (z ^ 0.8), and to study 
low mass clusters of galaxies (M 10 14 /i _1 M©). 



7 CONCLUSION 

We have presented the first shapelet analysis of weak grav- 
itational lensing surveys. We have constructed convergence 
maps of the CFHTLS Deep Dl field, and of 4 deg 2 of the 
CFTHLS Wide Wl field, which include the Dl field. We 
have detected six clusters of galaxies, through the lensing 
signal they ge nerate. Our Dl map is in good agreement with 
that of lGS07l . precedently created using the KSB shear mea- 



surement method. We combined our weak lensing data with 
the X-ray analysis of X MM-LSS CI cluster s lying in the 
same region of the sky l|Pacaud e t al. 2007). These three 
clusters catalogue are consistent. All our shapelet detec- 
tions have either an X-ray counterpart or a KSB detection. 
Counting our detections and accounting for the weak lensing 
selection function allowed us to constrain the power spec- 
trum normalisation <7 8 (fW0.24) ' 6 = 0.92^;?^. The com- 
bination of lensing masses and X-ray temperatures provided 
us with a new measurement of the mass-temperature rela- 
tion normalisation T» (or equivalently M*) for clusters of 
galaxies, M, = 2.7lt° 7 b \ lO 14 h _1 Af . Our results, though 
limited by low statistics and sample variance, are consistent 
with other current estimates. We also measured the slope of 
the mass-temperature relation, and found it consistent with 
self-similarity for low mass clusters, a = 1.60 ± 0.44. We 
have shown that one must measure both as and T* from 
combined weak lensing and X-ray surveys to investigate the 
discrepancy between independent measurements of as from 
different probes. 

Weak lensing surveys are becoming more and more ef- 
fective, and are currently being optimised for best extracting 
cosmological information. Optimal surveys will allow us to 
provide more accurate estim ates of as and T* , and to d isen- 
tangle the current as issue (|Amara fe Refregierj 12007m . We 
have compared the merits of weak lensing deep and wide 
blind surveys, based on the CFHTLS, at estimating <rg. We 
also looked at their merits at estimating T» while combined 
with an X-ray survey on their region of the sky. We found 
that for experiments driven by exposure time constraints, 
a wide survey will give ~ 3 times lower errors on the es- 
timates of both as and T*. To secure the measurement of 
as and M* with the current statistical accuracy, a 200 deg 2 
and a 300 deg 2 wide surveys will be needed respectively. 
We finally found that a 7000 deg 2 wide survey will be able 
to reach the 1% accuracy both on the power spectrum and 
mass-temperature relation normalisations. 
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